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ABSTRACT 

A Lagrangian, particle-based numerical method (tree code gravity plus smoothed par- 
ticle hydrodynamics) was used to simulate clump-clump collisions occurring within 
GMCs. The collisions examined were between TSM.Q clumps at a relative Mach number 
of = 3. The clumps were modelled using isothermal spheres which were individually 
in stable equilibrium. 

The collisions formed shock-compressed layers, out of which condensed approxi- 
mately co-planar protostellar discs of 7-6OM0 mass and 500-lOOOAU radius. Binary 
and multiple systems were the usual final state. Lower mass objects were also produced, 
but commonly underwent disruption or merger. Such objects occasionally survived by 
being ejected via a three-body slingshot event resulting from an encounter with a 
binary system. 

The impact parameter h denotes how offset the clumps are from one another, with 
low values corresponding to near-head on collisions, and high values corresponding to 
grazing collisions. Varying h altered the processes by which the protostellar systems 
formed. At low b a single central disc formed initially, and was then spun-up by an 
accretion flow, causing it to produce secondaries via rotational instabilities. At mid h 
the shocked layer which formed initially broke up into fragments, and discs were then 
formed via fragment merger. At large h single objects formed within the compressed 
leading edge of each clump. These became unbound from each other as b was increased 
further. 

The effect of changing numerical factors was examined by : (i) colliding clumps 
which had been re-oriented before the collision (thus altering the initial particle noise), 
and (ii) by quadrupling the number of particles in each clump (thus increasing the 
resolution of the simulation) . Both changes were found to affect the small-scale details 
of a collision, but leave the large scale morphology largely unaltered. 

It was concluded that clump-clump collisions provide a natural mechanism by 
which multiple protostellar systems may form. 

Key words: stars: formation - binaries: general - interstellar medium: clouds - meth- 
ods: numerical 



1 INTRODUCTION. 

This paper proceeds from the assumption that most star 
formation is dynamically triggered and produces clusters, 
rather than occuring spontaneously in a more widely dis- 
tributed mode. We explore the mechanics of star formation 
triggered by clump-clump collisions within molecular clouds. 
Such collisions can occur either within a single molecular 
cloud, as a result of the random velocities of the constituent 
clumps, or when two molecular clouds collide and the colli- 
sion is broken up into an ensemble of approximately simul- 
taneous clump-clump collisions. 

If the sub-structure within molecular clouds is hierar- 



chical, as has been inferred (Scalo 1985), the most common 
collisions will occur between clumps at similar levels of the 
hierarchy (and hence between c lumps of similar mass) . Ac- 
cording to perturbation theory (Whitworth et al. 1994), the 
characteristic masses of fragments condensing out of a shock- 
compressed layer are determined mainly by the post-shock 
temperature Ts- For typical values Ts ~ 10 to 30 K, the 
masses fall in the range 8 to 30 Mq. Therefore it is of par- 
ticular interest to consider collisions between clumps each a 
little more massive than this, since collisions involving much 
smaller clumps cannot produce fragmenting layers, whilst 
collisions between significantly larger clumps can probably 
be viewed as multiple collisions between smaller sub-clumps. 
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In this paper we have therefore restricted the parameter 
space to focus on collisions between equal mass clumps, each 
of mas s ~ 75M0 , and we have explored how the mechanics 



of fragmentation and star formation changes as a function 



of the impact parameter of the collision. 

The importance of magnetic fields in star formation is 
unclear. Given the low overall efficiency of star formation, 
it is possible that clump/clump collisions only trigger effi- 
cient st ar formation when the colliding clumps are threaded 



on the same flux tube. Under this circumstance, the critical 



surface density for gravitational condensation is not ampli- 
fied by the collision, and once the density is high enough the 
gas effectively decouples from the field, so the magnetic field 
has little influence on the events we have simulated. Alter- 
natively, the magnetic fleld may be destroyed by magnetic 
r econnection in certain clumps, in the manner described by 
( [Lubow fc Pringle 199^ ). 

§^ describes the simulation method and initial condi- 
tions, §^ discusses the results of collisions with varying im- 
pact parameter, §^ investigates the effect of numerical noise 
(which is introduced by the discretisation of the gas into 
particles) , §^ investigates the effect of increasing the resolu- 
tion of the simulation, and the overall conclusions are given 
in il 



2 METHOD. 

2.1 Numerical Details. 

Hydrodynamic forces were calculate d usin g smoothed parti- 
cle hvd rodvn amics ('SPH', see Lucv (1977), Gingold & Mon- 
aghan ( 1977 ), and Monaghan ( 1992 )). SPH is a Lagrangian 
method which represents a continuum fluid using a set of 
discrete sampling points. It is widely used in astrophysical 
problems because of its ability to handle large density con- 
trasts and complicated non-symmetric geometries. 

The sampling points are often referred to as 'particles' 
since they are assigned physical quantities and are moved ac- 
cording to the equations of motion of the fluid. A continuum 
quantity at arbitrary position, A{r), is found by performing 
a weighted sum using the value of A at each particle, with 
the weights depending on the distances from r. In fact, the 
weighting or 'kernel' function has flnite extent (the smooth- 
ing length 2h) so only neighbouring particles contribute. We 
allow the smoothing lengths to vary freely, both in time and 
from particle to particle, and choose to set them such that 
each particle sees, typically, 45- 55 neighbouring particles. 
We use the M4 weighting kernel ( Monaghan 1985 ), and the 
standard SPH artificial viscosity (Monaghan 198£) with pa- 



rameters a — 1 and f3 = 1 (this viscosity is necessary to 
prevent unphysical particle penetration in colliding fiows). 

Gravitational forces we re calculated using a Barnes-Hut 
tree (Barnes & Hut 1986), one of the tree-code methods 



which exploit the fact that at sufficiently large distances (de- 
termined by some tolerance parameter) a multipole expan- 
sion can provide a good approximation to the gravitational 
field of a group of particles. The force calculation is reduced 
from an 0{N^) process to one of 0{N log N). We choose to 
expand only up to the monopole term (effectively evaluat- 
ing the centre of gravity of gro ups of particles), but u se a 
tolerance parameter of 0=0.57 ( ^almon fc Warren 1994 ). 



To avoid spurious two-body events, gravitational inter- 
actions between particles at small separations mu st be soft- 
ened. Th e kernel-softening method was employed ( Turner et 
al. 1995). With this method the density profile of a particle 



is described by the kernel function (here, the M4 kernel), and 
interactions are unaltered outside a given softening length 
2e. We chose to use a constant e, so that gravity was strictly 
conservative. 



A leapfrog integration scheme was chosen (Hockney & 



Eastwood 1988), with the step-size allowed to vary from par- 
ticle to particle. However, these step-sizes were discretised 
by choosing a maximum size and then successively halv- 
ing; this allows the system of particles to be synchro nised . 
A similar approach was used by Navarro & White (199J) 
(though they employed a Runge-Kutta integration scheme). 
Since it is necessary to know the positions and velocities 
of all particles when performing an acceleration calculation, 
these quantities were linearly interpolated as required; in 
this manner the positions and velocities of all particles are 
known at multiples of the smallest step size in use. 

To ensure the code performed correctly, the following 
test simulations were successfully run : the free-fall collapse 
of a sphere of gravitating matter (a dynamic test of gravity), 
a 3D shocktube (a dynamic test of hydrodynamics) and the 
evolution of equilibrium isothermal spheres (an equilibrium 
test of both gravity and hydrodynamics). These tests are 



similar to those described in Turner et al. (1995) (it should 
be noted that the code used here was written in C and was 
developed independently from the Fortran code described in 
Turner et al). In addition, the self-similar evolution of an adi- 
abatic cloud of gas, as described by Evrard (1988), was also 



successfully modelled. This provides a dynamic test of both 
gravity and hydrodynamics and involves a variety of effects 
including rarefaction waves, freefall velocity flelds and stag- 
nation surfaces. By using appropriate boundary conditions, 
the similarity solution was reproduc ed m ore accuratel y and 
for much longer than either Evrard ( 1986 ) or Theuns ( 1992 ) 
achieved. 



2.2 Equation of State. 

We use a barotropic equation of state, i.e. the pressure is a 
function of density only. If we define the isothermal sound 

speed fliao by 



2 pkT 
m 



(1) 



(where for molecular gas m ~ 4 x 10 ^*g) then a^^^ is given 
by. 



+ 



((ag- a?) a?) (^)' 



P < Po; 



P> po; 



(2) 



where ao 



0.6km s ^, ai ~ 0.2km s ^, po — 1.2 x 



10"^°g cm-^ pi = 1 X 10""g cm-^ and 7 = 7/5. This 
equation of state may be described in terms of three regimes. 
Firstly, at low densities the gas is isothermal at ~100K. Sec- 
ondly, at intermediate densities the gas cools with increasing 
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density (T cx approximately), until T falls to -^lOK, 
and then the gas is approximately isothermal at ~10K. Fi- 
nally, at high densities the gas becomes adiabatic (because 
it becomes optically thick to its own cooling radiation) and 
r (X A smooth join is obtained between the last two 

regimes by taking their geometric mean. We believe that 
this equation of state mimics the broad properties of pro- 
tostellar gas which is optically thin at low densities (cool- 
ing radiatively with compression) , before becoming optically 
thick (whereupon it heats up adiabatically) . It should be 
emphasised that the adiabatic heating only becomes signifi- 
cant when the density has increased from its initial value of 



10"^^g cm-3 to > 10"" 



g cm 



2.3 Initial Conditions. 

The initial conditions are intended to represent two indepen- 
dent clumps within a GMC, just before a collision between 
the clumps takes place. Details of the complicated condi- 
tions within the GMC are omitted in favour of an idealised 
representation, which although simplified, captures the es- 
sential aspects of the clumps and the collision process. The 
clumps used are individually in stable equilibrium prior to 
the collision (with no imposed perturbations), and so any 
structure which subsequently forms must result purely from 
the collision, the equation of state of the gas, and any low 
amplitude seed noise introduced by the numerical method. 

The clumps are represented by spherically symmetric 
distributions of isothermal ideal gas in hydrostatic equilib- 
rium. Equating the gravitational and pressure forces for a 
symmetric distribution of gas gives 



Idp 
p dr 



GM{r) 



(3) 



Here M{r) = jATvr^p{r)dr (i.e. it is the mass interior to 
radius r). Using this, and writing the equation of state (|l|) so 
that dp/dr = af^^p d{\np)/dr, enables (^) to be re-arranged 
to the form 



dtp 

'di 



where the following substitutions have been made : 
i4nGp.)ir^ /p^\ 

fliso \ P J 



(4) 



(5) 



{po is the central density of the distribution). The density 
and mass profiles of such spheres may no w be obtained by 
numerical integration (see Chandresekhar ( 1939| ) and Turner 
et al (1995)). The resulting profiles are infinite in extent, so 



to produce a finite clump we specify the total mass of the 
clump and the isothermal sound speed, and we choose a 
value for ^i,, i.e. the value of ^ at the clump boundary, ^j, 
determines how centrally 'peaked' the density profile of the 
clump is. 

Because the clumps are truncated isothermal spheres, 
an external pressure must be exerted upon them to restore 
equilibrium. This is done by immersing them in an ini- 
tially uniform interclump medium (with density equal to 
that found at the edge of the clumps), and the interclump 
medium is in turn contained by a rigid boundary. To repre- 
sent this system, three types of particles are used. 'Clump' 



particles exert and feel forces due to both gravity and hydro- 
dynamics, and are free to move. 'Interclump' particles exert 
and feel hydrodynamic forces only, since they are intended to 
represent the rarefied regions which surround GMC clumps. 
These particles are also free to move. 'Boundary' particles 
exert hydrodynamic force but remain static. All particles 
have equal masses. The interclump region is chosen to be 
sufficiently large so that the results are not corrupted by 
edge effects propogating inwards. 

All these regions are produced from uniform, cubic par- 
ticle distributions. Such a distribution is obtained by placing 
particles at random positions within a cubic volume, apply- 
ing periodic boundary conditions, and allowing the system 
to evolve according to the hydrodynamic forces present. The 
evolution is terminated once the system has settled, i.e. the 
kinetic energy has decayed so that the distribution is static, 
and the density profile is uniform (within small scale dis- 
cretisation effects). The resulting particle positions are ap- 
proximately equidistant from each other, but not exactly 
regular (as on a grid or lattice). 

The interclump and boundary regions are then pro- 
duced from the uniform distribution by scaling it to obtain 
the required density, and joining and truncating copies of 
the resultant cubes to fill the required volume. 

A clump is constructed by extracting a sphere from the 
original uniform distribution and (non-linearly) scaling it to 
produce the required equilibrium density profile. It is then 
also allowed to settle again, to remove any small-scale un- 
balanced forces and motions resulting from the discrete rep- 
resentation of the gas. This is done by allowing it to evolve 
according to both hydrodynamic and gravitational forces for 
~ 20tff, within a spherical jacket of interclump and bound- 
ary particles. 

The final distribution is produced by removing two 
spherical volumes from the simulation interclump distribu- 
tion, and inserting both the settled clump, and a copy. The 
resulting initial conditions are as 'quiet' as possible, and free 
from grids or any other regular particle alignments, and any 
seed noise is of very low amplitude. 

Fig. ^ depicts the initial conditions schematically and 
shows the co-ordinate system used. The clumps are given 
velocities in the x direction, are aligned in the z direction, 
and may be offset in the y direction. This offset is described 
by the impact parameter b which is the vertical distance 
between the centres of the clumps ('b' in the figure) divided 
by the clump diameter. The right hand clump (at large x) 
is always higher (at larger y) than the left hand clump if 6 is 
non-zero. Hence b — 0.5 requires the lower edge of the right 
hand clump to be at the same height as the centre of the 
left hand clump. The distance between the centres of the 
clumps ('d' in the figure) is always set to be such that the 
two clumps are initially just touching. 

2.4 Simulations conducted. 

All collisions were between identical clumps with tempera- 
ture T — WOK (equating to a sound speed of a ~ 0.6km 
s~^), mass M = 75M0, and with ^f, = 3. This results in 
a clump radius of r ~ 0.6pc, and central, edge and mean 

densities Pcentrat = 1-1 X 10"^°g Cm~^, Pedge = 3.8 X 10"^^g 

cm~'^and p — 5.6 x 10~^^g cm~^. The freefall time of such 
a clump is tn ~ 0.9Myr. Each collision also took place at a 
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boundary 



interclump 




Figure 1. Clump-clump collision initial conditions. The clumps are immersed in an interclump medium which 
extends 4h beyond the bounding box containing the clumps, and this is in turn contained within a boundary layer 
2.2h thick. The clumps are offset by a distance 'b' in the y direction, and have a separation between centre-points 
of 'd'. 'd' is chosen such that the clumps are just touching. The clumps have equal, but opposite velocities in the x 
direction. 



relative velocity of ~ 1.8km s~^or a relative Mach number 
of = 3 (i.e. each clump was given a velocity of ~0.9km 
s~^), and each collision was followed for a period of IMyr. 

The effect of varying the impact parameter h was in- 
vestigated by a series of collisions spanning ti=0.1 to 0.7, in 
steps of 0.1. These collisions used 2000 particles per clump, 
and a gravitational softening parameter e=Q.002 pc. 

To investigate the effect of noise introduced by repre- 
senting the gas by a discrete set of particles, two additional 
collisions were performed at fe=0.5, but with the clumps re- 
oriented to change the particle positions within them. All 
other parameters remained as before. 

To investigate the effect of numerical resolution, a sin- 
gle high-resolution collision was also conducted with impact 
parameter 6=0.5, but employing 8000 particles per clump, 
and with a gravitational softening length of e=0. 00127 pc. 



2.5 Presentation of the results. 

The figures showing the results of the collisions are plots 
of column-density in which the grey-scale is logarithmic. In 
addition, contours which are equally separated in log-space 
are overlayed on these plots. 

It should be noted that the contours highlight relatively 
low density features. The compact condensations that form 
have much higher densities, but further contouring would 
be illegible. The high density of these protostellar conden- 
sations can be seen by examining the captions of the figures 
- these state the column-density range present. 

The results of the collisions are summarised in tables 
which also give the formation time of the first compact con- 
densation, and the values of various physical properties for 
the condensations formed. These properties are the radius 
(r*), the maximum density (n,), the mass (M,), and the 
magnitude of the specific angular momentum about the cen- 
tre of mass of the condensation (j*). The (non-specific) mag- 
nitude of angular momentum of a condensation (|J*|), and 
the ratio of rotational energy to self gravitational energy 
(/9,) can also be calculated. 

To determine these values, the constituent particles of 



the condensation must be found. This was done by discard- 
ing particles below a density threshold of p = 6.8 x 10~^^g 
cm~'^, and then associating each particle with its densest 
neighbour. The peaks of condensations may then be found 
by searching for particles which have themselves as densest 
neighbours, and the set of particles belonging to the conden- 
sation are those which point towards the peak directly, or 
indirectly through other particles. The values of the proper- 
ties listed above may then be calculated easily. 

The centres of mass may also be calculated, and per- 
forming this process over a time sequence of particle data 
produces the paths taken by the condensations. Figures of 
this type are used to illustrate the paths of objects produced 
after fragmentation. 

Since the collisions take place with clumps offset in 
the y-direction, and with velocities in the s-direction, the 
only non-zero component of angular momentum is the z- 
component, Lz- Since should be conserved, it can be used 
as an indicator of the accuracy of the code - the change in 
Lz for the entire system (i.e. summed across all particles 
within both clumps) should be small. The maximum per- 
centage change for each simulation is given in the summary 
tables. 



3 VARYING THE IMPACT PARAMETER. 

As the impact parameter b is altered, the mechanism by 
which protostellar discs are produced is changed. We con- 
centrate on examining results from the b=0.3, 0.5, and 0.7 
collisions. Results from all the collisions are summarised in 
Table |l|, and discussed in §3.4. 



The clumps were each represented by isothermal 
spheres of 2000 particles. The total number of particles used 
varied from 33755 for the 6=0.1 collision (2000 per clump, 
16519 representing the interclump, and 13236 forming the 
boundary) to 38896 for the 6=0.7 collision (2000 per clump, 
20284 representing the interclump, and 14612 forming the 
boundary). 
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Table 1. Summary of the standard collisions (see §|^). 

b t, r, n, M, j* AL^ 

(/dia) (Myr) (AU) (10"cm-3) (M©) (lO^^cm^s-i) (%) 



0.1 0.63 1330 16 60 2.1 -2 

One central object. 

0.2 0.62 1030 12 60 2.2 -2 

Six objects (2-6OM0), one dominant plus two smaller objects formed by rotational instability. 

0.3 0.64 880,1130 3.0,5.3 21,40 1.3,2.0 -2 

Three objects (3-40Mq), two dominant. Secondary formed by rotational instability from primary. 

0.4 0.68 580,860 4.0,2.7 13,27 0.5,1.5 -1 

Six objects (2-27Mq), two dominant. Secondary formed by lumpy accretion flow onto primary. 

0.5 0.78 890,920 1.2,1.5 17,17 1.3,1.3 -1 

Two objects, each formed by multiple fragment mergers, falling towards each other. 

0.6 0.89 840,830 0.5,0.8 11,11 1.2,1.1 -1 

Two objects, well separated. 

0.7 0.94 680,650 0.4,0.4 7,7 0.8,0.8 -1 

Two objects, well separated. 

r», M* and j» are determined using particles above a density threshold of n = 1.7 X 10** cm~^. t» is the earli- 
est time 2M0 of material forms a compact condensation, n, is the maximum density of a protostar at the end of 
the simulation. 



3.1 6=0.3 

The evolution of the system was dominated by accretion 
onto a central condensation and its subsequent rotational 
instability (Fig. ||). 

A single rotating spherical condensation rapidly formed 
at the collision centre, within a shock which was thin in 
the x-y plane, and somewhat extended in the z direction. 
The shock subsequently became more spindle-shaped (less 
extended in the z direction), and the condensation more 
disc-like. The smooth accretion flow along the shock caused 
the condensation to grow in size and mass, and increased its 
speciflc angular momentum causing it to became increas- 
ingly rotationally unstable. Material was then thrown off 
by the disc via weak spiral arms, and a surrounding difi'use 
envelope was produced out of which a distinct secondary 
condensed. This subsequently orbited the primary and grew 
in mass by 'mopping-up' the surrounding difi'use material, 
and by periodically intercepting the accretion flow. The fl- 
nal state consisted of a binary system of two protostellar 
discs, with a primary of 4OM0 and 1100 AU radius, and a 
secondary of 2IM0 and 880 AU radius. The final separation 
was 8820 AU, and the secondary completed two orbits of 
the primary. 

The evolution of various physical properties of the con- 
densations is shown in Fig. and is discussed below. The 
subscripts '1' and '2' will be used here to denote properties 
for the primary and secondary condensations respectively. 

Initially the accretion fiow within the shock-bounded 
spindle is focussed onto the central primary object, but a 
small ofiset rapidly develops, causing material streaming in 
from the upper-left to meet the lower edge of the conden- 
sation, and material from the lower-right to meet the upper 
edge of the condensation (Fig. ^ . The finite impact param- 




X / pc 

Figure 3. f)=0.3 collision. Offset accretion flow onto the central 
condensation at f=0.79 Myr. Crosses mark particle positions and 
lines are velocity vectors. Vectors are only plotted for particles in 
regions below a threshold density (~ 1 X 10~^^g cm""*), leaving 
the central condensation particle positions visible. 



eter of the collision causes the spindle to rotate about the 
angular momentum axis of the collision (the z-axis), and 
the accretion flow within the spindle evolves further and 
further away from the initial linear conflguration, becom- 
ing progressively more s-shaped and offset from the collision 
centre. Thus the angular momentum introduced to the cen- 
tral condensation is steadily increased, and the condensation 
is spun- up (|Ji| and Mi increase rapidly), causing it to be- 
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0.54- Myr. . 1 ' ' ' 1.01 iVtyr. 




Figure 2. 6=0.3 collision evolution. Accretion onto a central condensation is followed by rotational instabili- 
ties. [x=(-0.08,0.08)pc, j/=(-0.06,0.06)pc, contours at log^ (Af/cm-2) = (22.4, 22.7, 23.1, 23.4, 23.7), min=22.4, 
max=27.2.] 



come unstable increases). This can be seen in Fig. |^ at 
t <0.76 Myr. 

Once Pi - 0.25 and |Ji| ~ 1.4 x 10^^ g cm^ s'^ , ma- 
terial is shed via rotational instabilities (at t=0.76 Myr). 
Shedding and re-accretion cause oscillations in Mi and |Ji| 
until two fragments are permanently shed (f~0.77 Myr). It 
is worth noting that the instability does not just produce 
a simple spiral-arm structure. The process is more irregu- 
lar, as weak spiral arms are repeatedly formed and merge 
with diffuse material ejected earlier, and with material from 
the accretion flow. The two fragments produced are visible 
as the two traces below the label 'C in the mass evolution 
plot of Fig. ^. These subsequently merge to produce the sec- 
ondary (at t=0.86 Myr, the beginning of trace 'B' in Fig. 
%■ 

The increasing offset of the accretion flow causes in- 
coming material to swirl in with a large angular momentum 



relative to the discs. As the secondary orbits the primary it 
intercepts the flow and gains mass and angular momentum 
from it, so that M2 and | J2 1 increase, in addition some of the 
diffuse material around the primary is also collected. Mean- 
while, the primary is both accreting and occasionally shed- 
ding material from the edge of the disc (thus removing an- 
gular momentum, e.g. objects 'D' and 'E' in the mass evolu- 
tion graph in Fig. Mi increases slowly whilst |Ji| remains 
approximately constant. Thus f3\ (and hence the amount of 
rotational support) is decreased and the disc shrinks slightly 
(the 90% mass radius falls from 700 AU to 500 AU), whilst 
the peak density rises (continuing an almost linear trend 
since the formation of the primary). 

Fig. ^ demonstrates that the protostellar objects formed 
are discs, and that they are co-planar. The detailed internal 
structure of the discs is poorly modelled due to gravity soft- 
ening, and excessive shear viscosity, and so is not examined 
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Figure 4. b=0.3 collision. Evolution of condensation physical properties : mass M,, 100% and 90% radii r,, angular 
momentum about the condensation centre of mass |J»|, and the ratio of rotational to self gravitational energy /3». 
A minimum mass threshold of IOMq was applied to all except the mass plot. All graphs show the evolution of 
the primary (A), and the secondary (B). The mass evolution graph also shows the two fragments which merge to 
produce the secondary (C), a transient fragment thrown off and subsequently disrupted by the primary (D) and 
another fragment produced by the primary (E). 



(our main aim here is with the dynamics which form the 
discs rather than the subsequent evolution). 



3.2 6 = 0.5. 

The evolution of the system was dominated by fragmenta- 
tion and subsequent merger within the shock (Fig. P). Ini- 
tially, a single shock filament extending diagonally in the x-y 
plane was produced at the collision interface of the clumps. 
At about t=0.79 Myr, four fragments appeared in the shock 
in two close pairs, one above and to the left of the cen- 
tre, and the other below and to the right. Each pair subse- 
quently merged and accreted any further fragments which 
developed to produce a final state of two well defined proto- 
stellar condensations of similar size (~800 AU radius) and 
mass (I7M0), with a separation of 0.23 pc. These discs were 
falling towards each other along the shock at the end of the 



simulation, and would probably capture each other at peri- 
astron. 

Figure |^ illustrates the paths and mass evolutions of the 
fragments throughout the simulation, with corresponding 
positions and masses labelled. Fragments within the upper- 
left region are shown in the left-hand plots, and fragments 
within the lower-right region are shown in the right-hand 
plots. The evolution of both regions is broadly similar, with 
fragmentation initially producing objects of IM0 and lower 
(though these are equivalent to about 20-30 particles, and 
so are poorly sampled). Merger events are clearly visible, 
for example in the lower-right region fragments A2 and B2 
merge to produce C2 at 0.82 Myr, D2 forms at about the 
same time and merges with C2 to produce E2 at ~0.91 Myr, 
and finally fragment F2 appears and merges with E2 (at 0.97 
Myr) to produce the final disc G2 . These objects are also vis- 
ible at the right-hand side of the column density plots of Fig. 

H 

The paths taken by the fragments were affected by the 
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Figure 5. 6=0.3 collision. Protostellar discs at the end of the simulation (particle positions). Co-ordinates are 
centred on the centre of mass of the primary. Primary and secondary (left), and primary only (right). 



shock morphology. The shock initially formed at the colli- 
sion centre, and then propagated outwards as more of the 
clumps came into contact. Thus, the sequence of new frag- 
ments, i.e. those not produced by merger, similarly propa- 
gated outwards. At the lower-right region this sequence of 
fragments is A2, B2, D2, and F2. 



3.3 6 = 0.7. 

The collision produced a diagonal density enhancement (Fig. 

with the enhancement strongest at the ends of the di- 
agonal. These regions subsequently collapsed to produce 
two well separated condensations, and whilst the conden- 
sations were forming the leading edge of each clump ex- 
panded into the rarefied region behind the trailing edge of 
the other clump. When the simulation was terminated at 1 
Myr two condensations of mass 7M0 and radius ~660 AU 
were present, at a separation of 0.94 pc. This collision is in a 
relatively early stage of evolution, as the condensations form 
at later times than those in low impact parameter collisions. 



3.4 Dependence on the impact parameter. 

The main points of this discussion are summarised in Fig. 
% 

Low b collisions (6=0.2, 0.3) produce multiple objects 
by a process of accretion induced rotational fragmentation 
acting on an initi ally fo rmed central condensation (see also 
Whitworth et al (1995)). Unequal mass components result 



from this process, at least up to the end of the simulations 
presented here (1 Myr after the collision). Decreasing the im- 
pact parameter produces a higher mass primary (Mi=6OM0 
if 6=0.2, and 4OM0 if 6=0.3) since the reduced clump offset 



results in incoming material having less angular momentum. 
This slows the spin-up and allows more mass to accrete onto 
the primary. 

A collision of 6=0.4 includes mechanisms present both 
in the low 6 and 6=0.5 collisions, with the accretion flow 
onto the central disc eventually fragmenting. 

Higher 6 collisions (6=0.5, 0.6, 0.7) produce two equal 
mass protostellar objects in separate regions. The formation 
times and separations increase as 6 is increased. The systems 
are not well evolved by the end of the simulations, and the 
two objects in each collision have not interacted. However, 
the final objects in the 6=0.5 collision are bound, those in the 
6=0.6 case are marginally bound, and those in the 6=0.7 case 
are unbound. The objects in the two high impact parameter 
collisions will probably go on to accrete significant amounts 
of mass, and thus they may subsequently become bound. 
The ratio of the two condensation's mutual kinetic energy 
to gravitational potential energy is ~ 4.6 at the end of the 
simulation. 

The 6=0.5 collision produces its protostellar discs by 
a process of shock-induced gravitational fragmentation, fol- 
lowed by capture and merger between fragments (see also 
Turner et al ( |l99E| )). 

The high impact parameter collisions (6=0.6, 0.7) pro- 
duce objects at the leading edges of the clumps due to 
the compression of the impact triggering collapse. The col- 
lapse occurs because pressure support is reduced in the com- 
pressed regions as the gas cools. 

If 6a„itch is defined as the impact parameter below which 
rotational instabilities are dominant, then 6a„itch ~0.45. Un- 
equal mass systems result if 6 < 6switch, and equal mass 
systems (i.e. Mi ~ M2) if 6 > 6a„itch. It is expected that the 
binary systems produced when 6 < 6switch will have smaller 
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Figure 6. 6=0.5 collision evolution. Fragmentation in the shock is followed by merger between fragments. [x={- 
0.20,0.20)pc, j/=(-0.15,0.15)pc, contours at logi(,(Af/cm-2) = (21.9, 22.3, 22.6, 22.9), min=21.6, max=26.4.] 



separations than those when b > ^switch, since although the 
binary components of the latter are bound at formation they 
initially have much larger separations. 



4 VARYING NUMERICAL NOISE. 

Two additional collisions were performed which were iden- 
tical to the standard 6=0.5 collision, but with each clump 
re-oriented before the collision. This allowed the dependence 
of the results on particle noise to be investigated (the noise 
is a result of the discrete nature of the particle distribu- 
tion) . Since the clumps used have been relaxed, there are no 
grid-aligned directions, or any other preferential axes. Both 
clumps were rotated about an axis passing through their 
centres in the z-direction by 120° for one of the collisions, 
and by 240° for the other. The value b=0.5 was chosen be- 



cause it was thought to be the impact parameter most sen- 
sitive to noise. 

The final states of the original b-scan collision plus the 
two new orientations are shown in Fig. |l^, and are sum- 
marised in Table |^. Multiple fragmentation and subsequent 
mergers occurred in each collision, usually until two pairs 
of objects remained at each end of the shock. Merger be- 
tween each of the objects in the pairs was a possibility. 
At the right-hand side of the 120° case a 8M0 secondary 
merged with a I2M0 primary after completing about one 
orbit. Material was ejected during the merger process, and 
a single 17Mq object was produced. At the left-hand side 
merger did not occur, with a 6M0 object thrown into an el- 
liptical orbit around a I5M0 primary. More than one orbit 
was completed. Merger also did not occur at the right-hand 
side of the 240° case, where a 7M0 secondary encountered a 
I4M0 primary. Instead, momentum was exchanged and the 
primary was ejected towards the collision centre. Merger of 
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Figure 7. b=0.5 collision. Evolution of condensation masses and positions (position markers scale linearly with 
mass). Condensations appearing to the left and above the centre of collision (left) are separated from those appearing 
to the right and below (right). Each fragment is labelled so that the corresponding masses and positions are apparent. 



Table 2. Summary of the re-oriented clump collisions (see §y). 

b t, r* n* M-t j, ALz 

(/dia) (Myr) (AU) (10"cm-3) (Mq) (lO^^cm^s-i) (%) 



0.5 0.78 890,920 1.2,1.5 17,17 1.3,1.3 -1 

Rotated by 0° (standard collision). Two objects, each formed by multiple fragment mergers. 

0.5 0.76 890,780 1.1,1.5 15,17 1.2,1.3 -2 

Rotated by 120°. Three objects (6-17Mq), two dominant formed by multiple fragment mergers. 

0.5 0.76 790,750 0.9,1.3 12,14 1.1,1.0 -2 

Rotated by 240°. Four objects (8-14Mq), formed by multiple fragment mergers. 

r«, Mt and j* are determined using particles above a density threshold of n = 1.7 X 10* cm~^. is the earli- 
est time 2M0 of material forms a compact condensation, n, is the maximum density of a protostar at the end of 
the simulation. 



small fragments occured at the left-hand side to produce 
a I2M0 object, and also just left of the centre producing 
another 8M0 object. 

The shock size, shock shape, time of fragmentation, the 
total mass of discs produced, and the regions in which these 
discs appeared were broadly similar across all the collisions. 
The differences concerned where fragmentation occurred, 
and whether final merger events took place. Fragmentation 
was usually localised within two regions in each collision. 
These areas were small enough to allow the objects pro- 



duced within each region to interact, and to merge until the 
last pair of objects remained. By this stage the objects had 
evolved into discs, and the probability of merger was found 
to be ~50%. In one of the collisions however, fragmentation 
also occurred in a more central region, and a single object 
was produced. The central disc was at a sufficient distance 
from the other fragmentation regions to preclude interaction 
up to IMyr after the time of collision. 
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Figure 8. b=0.7 collision evolution (upper left region shown, lower right is very similar). The collision perturbs 
and compresses the leading edge of each clump, causing collapse. [a;=(-0.60,0.00)pc, y=(0.00,0.45)pc, contours at 
logiQ(Ar/cm-2) = (21. 6, 21.9, 22.3, 22.6), min=-oo, max=25.9.] 



5 VARYING NUMERICAL RESOLUTION. 

To investigate the effect of resolution on the outcome of a 
simulation, the standard 6=0.5 collision was repeated with 
the number of particles in each clump increased by a fac- 
tor of four. This increased the linear spatial resolution by 
a factor of 43, or ~1.6. Accordingly, the gravity smoothing 
length e was dropped by the same factor, from 0.002 pc to 
0.00127 pc. A similar change was not required for the hy- 
drodynamic smoothing length h, since the size is controlled 
by the requirement that a reasonable number of neighbours 
be found within each kernel (this automatically reduces h to 
a suitable value). A total of 95411 particles were used (7999 
per clump, 52031 representing the interclump, and 27382 
forming the boundary). 

The collision evolution followed a similar pattern to the 
equivalent standard-resolution collision, but the shock pro- 



duced was thinner, and the fragmentation was better re- 
solved. The result is summarised in Table ^ 

Fragmentation initially occurred near to the collision 
centre within the shock bounded region. As the collision 
progressed, the density in regions of the shock further from 
the centre increased, and fragmentation spread out along the 
shock. By 0.92 Myr fragment merger had produced an ap- 
proximately symmetric configuration with two pairs of discs 
on either side of a small IMq central condensation. The in- 
ner discs formed from mergers close to the collision centre, 
and had masses of IIM© each. The outer discs formed from 
mergers in the outer regions of the shock, and were slightly 
smaller (each had a mass of 8M0). This configuration is sim- 
ilar to that found in the various standard-resolution 6=0.5 
collisions before the final encounters at either side of the 
shock take place. Here, the two small discs encountered the 
larger discs at 0.95 Myr and merger did not occur in either 
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Figure 9. Summary of the results and the disc formatio n rn echanisms which occur in the clump-clump collision 
simulations, as a function of the impact parameter b. See § [j.4| for a fuller discussion. 



Table 3. Summary of the high-resolution collision (see [^). 

b t, r* n* M* j, AL^ 

(/dia) (Myr) (AU) (10"cm-3) (Mq) (lO^^cm^s-i) (%) 



0.5 0.77 560,580 3.4,4.7 11,12 0.7,0.8 -2 

Six objects (I-I2M0). Two binaries (9,12 & 11,12M0), components form by mult. frag, mergers. 

r«, M» and j* are determined using particles above a density threshold of n = 1.7 X 10* cm~^. is the earli- 
est time 2M0 of material forms a compact condensation, n* is the maximum density of a protostar at the end of 
the simulation. 



encounter. Instead, matter was ejected via the production 
and subsequent detachment of tidally induced arms (which 
also removed angular momentum and increased the binding 
energy of each pair of discs), and binary systems resuhed. 
The final state of the system is shown in Fig. |l^ (left). 

When compared to the standard-resolution simulations, 
some differences of detail are apparent : the binary com- 
ponents are of slightly smaller mass (9-12M0 instead of 
I2-I7M0), and more significantly the objects have higher 
peak densities, less specific spin angular momentum, and 
are smaller (at <600 AU instead of 800-900 AU). Thus the 
objects produced in the high-resolution simulation are more 
tightly bound (a consequence of the reduction in e and h). 
This reduction in size and increase in density is apparent 
throughout the simulation and causes the discs in the final 
encounters to resist merger, and to favour capture. Since 
the high resolution simulation is more realistic, it is expected 
that the production of binaries would be favoured in nature. 

A plot of fragment paths reveals many more objects 
in the high-resolution simulation. However, since a density 
threshold is applied, this means that denser fragments are 
more common. If in addition a IM0 threshold is applied 
(Fig. right) then the plot closely resembles the standard- 
resolution plots (Fig. these are left almost unaltered by 
applying the mass threshold) . Thus increasing the resolution 
allows better simulation of early small and dense fragments, 
but these subsequently merge and the larger scale behaviour 
is little changed. 



Another small scale difference was the ability to model 
late forming low-mass fragments appearing at the outer 
edges of the shock, far from the collision centre. These frag- 
ments were few in number, and their subsequent encounters 
with the binaries usually resulted in disruption or merger 
with one of the binary discs. However, one of the low-mass 
fragments encountered the two left-hand discs just as they 
were at periastron during their first encounter. A three-body 
event occurred and the IM0 fragment was ejected. Fig. [13 
shows a time sequence of this encounter showing the ejec- 
tion of this small object, and the tidal arms excited in the 
main discs. This is a process by which low mass objects may 
survive, potentially creating lower mass stars. 



6 CONCLUSIONS. 

The parameter space investigated was restricted to clumps 
of equal mass. We believe that relaxing this requirement to 
consider collisions between clumps of similar mass (which 
probably encompasses the majority of collisions occurring 
in giant molecular clouds), will produce results not greatly 



different to those presented here. Larson's relations (Larson 



1981) suggest that clump radius and density vary slowly 
with mass (r oc M^^^, p oc M~^^^), and so the geometry 
of such collisions will be largely unaltered. However, the re- 
sult of any collisions that occur between clumps of widely 
differing mass may well be significantly different, since both 
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Figure 10. Re-oriented, b=0.5 collisions. Final states : original collision (top), clumps re-oriented by 120deg about 
the z-axis (middle), and by 240deg about the 2-axis (bottom). Left : column density. [i;=(-0.20,0.20)pc, y={- 
0.15,0.15)pc, contours at logio(Ar/cm-2) = (21.9, 22.3, 22.6, 22.9), min=21.6, max=26.4.] Right : paths taken by 
condensations throughout the simulation, and final objects shown as circles with area proportional to mass. 



0,2 



the geometry and the extent of the shocked region will be 
substantially changed. 

The internal properties of the discs produced here, such 
as density and velocity profiles, cannot be determined reli- 
ably. This is due both to the excessive shear viscosity intro- 



duced by the SPH artificial viscosity formulation, and the 
comparatively large size of e required to prevent infinite den- 
sity collapse. A small e also results in a prohibitively small 
timestep). 

Unfortunately, it is also not possible to follow the sub- 
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Figure 11. High-resolution, b 


=0.5 collision. Final state. Left 


column density. [x={-0.20,0.20)pc, ?/={-0.15,0.15)pc. 


contours at log]^g(Af/cm~^) = 


(21.9, 22.3, 22.6, 22.9), min=21.6, max=26.7.] Right : paths taken by condensations 



of IMq and above throughout the simulation, and final objects shown as circles with area proportional to mass. 
Although the high-resolution simulation captures more detail of early low-mass objects than the standard-resolution 
simulation, the paths taken by subsequent objects (shown here) are similar (though final mergers do not occur and 
binaries are the end result). 



sequent evolution of the protostellar objects formed in these 
simulations. This evolution seems likely to involve a compe- 
tition between the internal evolution of the individual pro- 
tostellar discs (which converts a massive protostellar disc 
into a central star plus an emaciated circumstellar disc) and 
impulsive interactions between the discs (which may cause 
them to fragment and spawn additional protostellar discs in 
closer binary systems). In order to explore this possibility a 
new SPH viscosity formulation has been developed in which 
the bul k and shear compon ents can be regulated indepen- 
dently (Watkins et al. 1996). Simulations which concentrate 
on examining solely the interaction between a protostellar 
disc and either a naked star or another protostellar disc, hav e 
been performed using this method (Watkins et al. 1998). 
Such interactions are found to be effective at spawning new 
protostars and closer binary systems. 

The larger scale processes which deliver the material 
that produces the discs in the simulations presented here are 
believed to be well resolved. Thus, it seems that the mech- 
anisms which produce the discs are viable, and processes 
such as fragmentation and accretion-induced rotational in- 
stabilities are well modelled. The relatively large value of e 
(compared to the typical disc radius) produces discs which 
are more susceptible to disruption than they would be other- 
wise (since they are less well bound). The presence of many 
surviving discs, in spite of the high e, indicates that disc for- 
mation is a likely consequence of a clump-clump collisions in 
nature. In addition, quantites such as the masses, positions 
and orientations of the discs are expected to be unaffected 
by numerical factors, and so are believed to be reliable. 

Low impact parameter collisions (fe=0.2, 0.3) produced 
binary or multiple systems by the rotational instability of 
an initially formed central disc. In these systems the com- 
ponents tended to have unequal masses. Mid-b collisions pro- 
duced two bound objects of similar mass which formed via 
repeated mergers of fragments which appeared in the shock- 
compressed layer produced at the collision interface. High-fe 



collisions (b=0.6, 0.7) produced discs, also of similar mass, 
at the leading edges of each clump - these were unbound in 
the 6=0.7 case but were still accreting significant amounts 
of matter at the simulation end (1 Myr after collision). If 
the simulations could be followed further, a larger fraction 
of the mass would end up in protostellar discs. However, this 
is because we have not included any feedback processes, of 
the sort which might cut off accretion onto a protostellar 
disc. 

The numerical noise introduced by discretising the gas 
into a system of particles did not appear significantly to af- 
fect the results. Collisions between re-oriented clumps pro- 
duced results broadly similar to the original collision simu- 
lation. 

An increase in resolution also left the results of a sim- 
ulation largely unchanged - early small scale fragmentation 
was better resolved but subsequent evolution was left largely 
unaltered, though the final discs were slightly smaller and 
of slightly higher density. 

Although not presented here, collisions between spin- 
ning clumps were also investigated. The effect of a spin rate 
similar to that observed for such clumps was found to be 
negligible - the kinetic energy and angular momentum as- 
sociated with the spins of the individual clumps are much 
smaller than the translational kinetic energy and the orbital 
angular momentum associated with the collision. 

It is concluded that clump-clump collisions provide a 
realistic mechanism by which protostars may be formed. In 
fact, binary or multiple systems are often the result; this 
is in good agreement with observations of star-forming re- 
gions, which similarly find a preponderance of multiples and 
binaries. Collisions with a wide range of impact parameters 
were found to produce two or more protostellar objects. 
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Figure 12. High-resolution, 6=0.5 collision. Encounter between two discs (8 and 11 M©), and a IMq fragment. The 
fragment is ejected from the three-body system to leave a binary system. [a;=(-0. 165,-0. 097)pc, j/=(0.006,0.074)pc, 
contours at Iogio(Ar/cm-2) = (22.2, 22.7, 23.2, 23.7).] 
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